arXiv: 1504.01290v2 [cond-mat.stat-mech] 22 Apr 2015 


Glass transitions and scaling laws within an 
alternative mode-coupling theory 

Wolfgang Gotze 1 and Rolf Schilling 2 

1 Physik Department, Technische Universitat Munchen, 

James-Franck-Str. 1, D- 85747 Garching, Germany 

2 Institut fur Physik, Johannes Gutenberg-Universitat, 

Staudinger Weg 9, D-55099 Mainz, Germany 
(Dated: April 23, 2015) 

Idealized glass transitions are discussed within a novel mode-coupling theory (TMCT) proposed 
by Tokuyama(Physica A395,31(2014)). This is done in order to identify common grounds with 
and differences to the conventional mode-coupling theory (MCT). It is proven that both theories 
imply the same scaling laws for the transition dynamics, which are characterized by two power-law 
decay functions and two diverging power-law time scales. However, the values for the corresponding 
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parts of the correlation functions. It is also demonstrated for a schematic model that the TMCT 
neither leads to the MCT scenarios for transition-line crossings nor for the appearance of higher- 
order glass-transition singularities. 
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I. INTRODUCTION 


With decreasing temperature T or increasing density 
p the dynamics of a liquid slows down drastically, and 
finally it arrests in a glass state. This holds provided 
crystallization can be avoided. An outstanding feature 
of glassy dynamics is stretching: the relaxation spectra 
extend over many orders of magnitude for the frequency 
variations. One cannot describe meaningfully the relax¬ 
ation functions as superpositions of a number of expo¬ 
nential decay functions. Mode-coupling theory (MCT) 
is an attempt to describe the indicated features within 
the microscopic theory of simple liquids [1, 2]. The ba¬ 
sic quantity of this approach is the intermediate scatter¬ 
ing function S(q,t) = (p(q,t)* p(q)}/N. For an N-particle 
system it describes the variations of the correlations with 
increasing time t of density fluctuations p{q, t ) of wave 
vector q ., \q\ = q. The equation of motion reads 


S{q, t)+v q S{q, t)+n 2 q [S(q, t)+ 


r 

o 


m q (t—t')S(q,t')dt'] = 0 . 


(1) 

It is to be solved for the initial conditions S(q, t = 0) = S q 
and S(q, t = 0) = 0. This formula generalizes a Navier- 
Stokes-equation result by replacing the sound frequency 
by D q and by complementing the friction term u q by a re¬ 
tardation term. The latter is quantified by a fluctuating- 
force kernel m q (t) [3, 4]. Here S q denotes the static struc¬ 
ture factor, which depends smoothly on T and p. Often it 
is more convenient to discuss dynamics in the domain of 
complex frequencies z, Im(z) > 0. Eq. (1) is equivalent 
to the double-fraction formula 


S(q, z)/S q = -l/[z - D 2 q /[z + iv q + n 2 m q (z)]\ . (2) 


Here and in the following Laplace transforms F(z) of 
functions F(t) are defined by 

n OO 

F(z) = LT[F{t)\(z) = i / F(t) exp{izt]dt . (3) 

Jo 

Eqs. (1) and (2) are still exact, but with an unknown 
kernel m q {t) and rn q (z), respectively. To close the equa¬ 
tion of motion, an ansatz for the kernel was motivated, 
which expresses m q [t) as a polynomial of S(q,t). De¬ 
tails will be specified in Sec. IIA. Analogous MCT equa¬ 
tions have been formulated for other functions, like the 
density-fluctuation correlator S^ s \q,t) for a tagged par¬ 
ticle moving in the liquid [1, 2]. 

MCT describes ideal glass transitions occurring at 
some critical temperature T c or critical density p c . For 
T > T c or p < p c the long-time limits of S(q,t) van¬ 
ish as expected for a liquid state. However, for T < 
T c or p > p c , the correlations arrest for long times: 
limt-^oo [S (q, t)/S q ] = lim^o ~[zS{q, z)/S q ) = f q . The 
positive arrested parts f q are the Debye-Waller factors of 
an amorphous solid, i.e. of a glass. 

MCT also describes the evolution of a stretched dy¬ 
namics as precursor of the ideal glass transition. Uni¬ 
versal formulas have been derived in the limit T tending 
to T c or p approaching p c . Two scaling laws describe 
stretching as interplay of two power-law-decay functions. 
The subtle dependence of the dynamics due to changes 
of the small parameters (T — T c ) or ( p c — p) is caused by 
two power-law-divergent time scales, which are specified 
by two anomalous exponents [2]. 

So far it has not been possible to control the short¬ 
comings of the MCT ansatz for the kernel m q (t). This 
has challenged extensions and modifications of the MCT 
[5-10] . Particularly stimulating is the work of Biroli and 
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Bouchaud [10], which is built on equations for averages of 
products of four density fluctuations. They conclude that 
MCT is a mean-field approximation of their more general 
equations. But studies of the MCT for hard-sphere sys¬ 
tems in high-dimensional spaces do not provide definite 
validity of this conclusion [ 11 , 12 ], 

Recently, an alternative to MCT has been studied by 
Tokuyama [13] which will be called TMCT. While MCT 
is imbedded in the Zwanzig-Mori projection-operator for¬ 
malism [3, 4], the alternative TMCT is based on a differ¬ 
ent projection-operator formalism [14]. Like in MCT also 
in TMCT an ansatz for a relaxation kernel m q (t ) is re¬ 
quired. It was suggested [13] to use the same expression 
as was motivated in MCT. 

It was shown already that TMCT overestimates the 
glass formation less than MCT does [15]. It was also 
proven that the non-Gaussian parameter for t = 0 van¬ 
ishes within TMCT while MCT yields a non-zero value 
[13] . The aim of the present paper is to investigate within 
TMCT the glass transition scenarios and the long-time 
dynamics and to compare with the corresponding results 
within MCT, thereby providing a motivation to test both 
theories by experiments and molecular-dynamics simula¬ 
tions. 

The present paper is organized as follows. Sec. II 
presents the equations of motion for MCT and TMCT, 
reviews the basic concepts and manipulations in order 
to discuss the glass-transition singularities and it con¬ 
tains the proof that TMCT - contrary to the MCT - 
does not exhibit continuous glass transitions. Sec. Ill 
and Sec. IV discuss within TMCT the validity of the first 
and second scaling law, respectively. The glass-transition 
diagram of a schematic model is calculated analytically 
within TMCT and its properties are compared with those 
calculated within MCT in Sec. V. The final section VI 
presents a summary and conclusions. In order to keep 
the presentation of our work self-contained the crucial 
concepts and steps will be explained within each section. 
Some technical details are put into an appendix. 


II. GLASS-TRANSITION SINGULARITIES 

In a first step, the equations of motion for MCT and 
TMCT will be presented in subsection A. Subsections B 
and C will describe necessary tools needed for the dis¬ 
cussion of the glass-transition singularities and the dy¬ 
namical features. Based on this, in subsection D it will 
be proven that TMCT can not describe continuous glass 
transitions. 


A. Equations of motion 

The dynamics shall be described by the normalized 
correlation functions cf> q (t ) = S(q,t)/S q . The wave num¬ 
bers shall be discretized to a set of M values for the 


moduli , - - - ,qM- It is the goal to evaluate the ar¬ 
ray 4>(t) = where k cor¬ 

responds to q k . This shorthand notation will be used 
throughout the following. As mentioned above, arrays of 
kernels m(i) have to be defined. This is done by intro¬ 
ducing arrays !F[P, x] of mode-coupling polynomials: 

no 

Wx]=£ £ v q { t.. kn (P)x kl ---x kn . (4) 

n— 1 k\---k n 

P = (Pi,--- , Pd) denotes a point in the d-dimensional 
manifold of control parameters. These points specify 
the equilibrium state of the system. The coefficients 
Vill. kn (P) play the role of the coupling constants. They 
are non-negative smooth functions of the control param¬ 
eters. The ansatz for the mode-coupling kernel reads 

m(t) = F[P, 4>(t)\ . (5) 

For studies of simple liquids quadratic polynomials have 
been motivated and the (P) are expressed by S q 
[1], For a study of a spin-glass-transition model only 
coefficients for n = 1 and n = 3 have been used [16]. 
Coefficients for n = 2 had to be ignored because of sym¬ 
metry reasons. In a study of the percolation problem for 
a hard sphere moving in a random arrangement of point 
scatterers only linear mode-coupling polynomials occur, 
i.e. I ^q' r kl ... kn (P) = 0 for n > 1 [17]. Equation (4) rep¬ 
resents a generalized mode-coupling polynomial with a 
lower and an upper cut-off for q at q\ and Qm-, respec¬ 
tively, and with discretized wave numbers. Since, e.g. 
for a liquid, the coefficients V q n kl ... kn {P) involve the di¬ 
rect correlation functions, they decay for the wavenum¬ 
bers going to infinity. This decay is such that for M and 
qM large enough the solutions of the MCT equations are 
practically independent of M and qM- The properties of 
a discontinuous glass transition are also robust under the 
continuum limit of the wave vectors, including the limit 
qi —> 0. However, for continuous transitions, e.g. for the 
tagged particle dynamics, the continuum limit is more 
subtle, since it also involves the limit qi —> 0. In Ref. [18] 
it was proven that in three spatial dimensions the transi¬ 
tion survives in that limit. But, the transition point and 
properties close to it, e.g. the dependence of the arrested 
parts of the corresponding correlators on the distance to 
the transition point, depend weakly on < 71 . In two dimen¬ 
sions the transition point, i.e. the critical density, shifts 
to zero for qi —> 0 . 

The MCT, as well as the TMCT, are based on two 
equations of motion. The first one is shared by both the¬ 
ories. It expresses a current correlator K(z) in terms of 
the Laplace transform of m (t) as Zwanzig-Mori fraction: 

K q (z) =-Vl 2 q /[z + w q +Vl 2 q rh q (z)\ . ( 6 ) 

The frequencies Q q and v q specify the transient dynamics. 
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The second equation of motion formulates a connec¬ 
tion between the density-fluctuation correlators and the 
current correlators. The MCT uses again a Zwanzig-Mori 
fraction: 

MCT : 4> q (z)=-l/[z + K q (z)\ . 

(7a) 

Combining Eqs. ( 6 ) and (7a) one gets Eq. (2). For the 
TMCT one appeals to an alternative projection-operator 
formalism [14] in order to get a connection between K(f) 
and in form of a first order differential equation 
[13]: <p q (t) + [Jq K q (t')dt'](j) q (t) = 0. The latter has to 
be solved for the initial condition <j> q (t = 0) = 1. One 
gets - In {(j) q {t)) = E q {t) := f*(t - t')K q (t')dt', i.e. 

TMCT : 0,(i) = exp(-£7,(i)); E q (z) = - K q {z)/z 2 . 

(7b) 

Two comments on the basic equations of both theories 
might be added. First, in order to use a unified descrip¬ 
tion for MCT and TMCT notations have been changed 
relative to the preceding literature on TMCT. This con¬ 
cerns the damping constant and the memory kernel. The 
microscopic frequency term(first term on the r.h.s. of 
Eq.(2) in Ref. [15]) has been denoted by H 2 , which is equal 
to K q {t = 0). The major point, however, is that the func¬ 
tion K c (q, t ) from Eq. (2) of Ref. [15] has been denoted by 
E q (t). Taking this into account, the Laplace transform 
of that equation is identical to our Eq. ( 6 ), and Eq.(l) 
in Ref. [15] coincides with Eq. (7b) of the present contri¬ 
bution. Second, since both theories are based on differ¬ 
ent projection formalisms, the current correlator K q (t) 
in MCT and TMCT are different functions. 


B. Fixed-point equations 

As mentioned in the introduction, glass states P are 
characterized by the arrest of density-fluctuation corre¬ 
lations: f = limt-xxj 4>(t) = lini z _ ) .o [—z<f>{z)\. Equa¬ 
tion (5) implies arrest of the mode-coupling kernels: 
T[P, f] = lim^oo m(i) = lim^ol - 2 m(z)]. Formula ( 6 ) 
leads to an equation for the low-frequency asymptote of 
the current correlators: 

lim[K q (z)/z\ = l/E q [P,f], q = 1, - ■ ■ ,M . ( 8 ) 

z—>0 

Substitution of this formula into Eq. (7a) yields f q = 
1/[1 + l/JF q [P, f]]. Using Eqs. (7b) and ( 8 ) one arrives at 
fq = exp [lim z -to[zE q (z)]\ = exp[-lim z _>o [K q {z)/z\] = 
exp[— l/F q [P, f]] . Defining mappings T q [P, x] of the man¬ 
ifold of arrays x into itself by 


the arrested parts f are fixed points of T q [P, x]: 


f = T[P, f] ■ (10) 

For the manifold of arrays a semi-ordering can be de¬ 
fined by: x > y if x q > y q , q = 1, • • • , M. Introduc¬ 
ing the special arrays 0 = ( 0 , • • • , 0 ) and 1 = ( 1 , • • • , 1 ) 
one can write 1 > f > 0. There hold the monotony 
properties: if x > y > 0 one gets P\P 1 x] > F[P,y\ 
and T[P 1 x] > T\P, y]. A sequence of arrays g^ n \ 
n = 0 , 1 , - - • can be defined recursively by: g(°) = 1 , 
gO-H) _ 7~[P, g(")]. Since g(") > g(" +1 ) > 0 , the se¬ 
quence converges towards some fixed point f* > 0. 


C. The stability matrix 

The analysis in this subsection follows that perfomed 
for MCT. For more details the reader may consult 
Secs. 4.3.3 and 4.3.4 of Ref. [2]. The fixed-point equa¬ 
tion Eq. (10) formulates M implicit relations T q [P, f] = 
f q — T q [P, f] = 0 , q = 1, ,M for the M compo¬ 
nents of a fixed point f (P). Let us assume that P* 
is a glass state specified by the fixed point f(P*) = 
f*. There appears the question whether this special 
solution can be extended to a neighborhood of states 
P “near” P* . The implicit-function theorem answers 
this question affirmatively provided the M-by-M ma¬ 
trix with elements T q ^[P*, f*) = dl q /dxk[P*, f*] can 
be inverted. Writing I q ^[P, x] = < 5 qi k — A q ^(P,x) with 
A q ^[P, x] = dT q [P, x]/dxk one concludes the follow¬ 
ing. If the matrix (A 9 i fc[P, x]) does not have an eigen¬ 
value unity, a smooth continuation f(P) of the fixed 
point array f* is defined uniquely for a whole neigh¬ 
borhood of states P near P*. The eigenvalue condi¬ 
tion for (Aq t fe[P, x]) can be replaced by one for matrices 
{A q ,k[P, x]) = (s g (x) _ 1 A g , fc [P,x]s fc (x)), s q (x) > 0, since 
both matrices are connected by a similarity transforma¬ 
tion. The choice of s(x) is a matter of convention. Here 
we will choose s q (x) such that s q (f(P))P q [P, f(P)] = 1. 
This choice simplifies the calculations in Sec. III. The 
new matrix (A qt k[P, x]) is called the stability matrix of 
the theory. 

From Eqs. (9a) and ( 10 ) one gets A q ^[P, x] = 
(1 — x g ) 2 9P g [P, x]/dxk and Eqs. (9b) and (10) yield 
A qt k[P, x] = x g (ln(l /x q )) 2 dP q [P,x\/dx k . Use of 
Eqs. (9a), (9b), (10) and the above convention for s q (x) 
leads to 


MCT : S q (x) = (1 - Xq) / X q , 
TMCT : s q (x) = ln(l/a; q ) . 

Then the stability matrices are given by 


(9a) 

(9b) 


MCT : T,[P,x] = 1/[1 + l/P 9 [P,x]] , 
TMCT: T q [P,x]=exp[-l/P q [P,x]] , 


( 11 a) 

(lib) 


4 


MCT : A gifc [P,x] = x q (l - x q )dP q [P,x]/dx k (l - x k )x k l , 

(12a) 

TMCT : A q} k[P,x] = x q \n(l/x q )dP q [P, x)/dx k ln(l/xfe) . 

(12b) 

Generically, the M 2 elements of the stability matrix are 
positive. Some consequences of this property shall be 
cited, which will be needed in the following. The matrix 
has a non-degenerate maximum eigenvalue to be denoted 
by e(P). All other eigenvalues have a modulus smaller 
than e(P). The implicit-function theorem permits the 
following conclusions. Let P* denote a state with a fixed 
point f* and e(P*) < 1. Then there exists a neighbor¬ 
hood of states P with fixed points f(P) and e(P) < 1. 
The M components of f(P) and e(P) are smooth func¬ 
tions of the control parameters and f(P*) = f*. Such P* 
are called regular glass states. 

Glass-transition singularities or critical states P c with 
arrested parts f c = f(P c ) are characterized by a critical 
stability matrix (A q k ) = (A q ^[P c ,f c ]) with maximum 
eigenvalue unity: e(P c ) = 1. Let a* and a denote the 
left and right eigenvector of (A qk ), respectively: 

a pAp k =a k > ^^g,p a p = a g- (13a) 

p p 

Generically all components of the eigenvectors can be 
chosen positive: 

a k > 0, a q > 0 ,k,q = l,---,M . (13b) 

The eigenvectors shall be fixed uniquely by the conven¬ 
tions 

a q a 1 = 1 > a q a <l a <l/fq = 1 ' ( 13c ) 

9 9 

In Secs. Ill and IV there will appear the problem to 
understand the solution X of a set of M linear equations. 
These are specified by the critical stability matrix and 
by an array I of inhomogeneities: 


XX* - Alk]X k = I q , q = (14a) 

k 

The solubility condition for this problem reads 

9 J 9=° ' ( 14b ) 

9 

The general solution of Eq. (14a) is the sum of two terms. 
The first one is an arbitrary multiple X of the components 
of a. The second one is a linear combination of the com¬ 
ponents of I with coefficients R q ,k given by the elements 
of the stability matrix: 

X, = Xa q + X q ■ X q = Y J Rqkh ■ (14c) 

k 


D. Absence of continuous glass transitions within 
TMCT 

Let P e , —eo < e < £o, denote a path in the manifold 
of control parameters. It shall be a smooth one in the 
sense that all mode-coupling coefficients I / ^".j... fcn (P e ) are 
smooth functions of the path parameter e. The conven¬ 
tion shall be made: states for e < 0 are liquids, the ones 
for e > 0 shall denote regular glass states, and P e= o = P c 
is a glass-transition singularity. It will be the goal of this 
subsection to prove the theorem: 

The TMCT equations (9b) and (10) with the mode¬ 
coupling polynomial (4) cannot yield arrested parts f(P e ) 
which change continuously with changes of e at P e —o = 
P c . 

The arrays x which enter the mode-coupling polyno¬ 
mial obey the restriction 1 > x > 0. There is only a 
finite number of terms x kl ---Xk n which enter 

Eq.(4). Hence there is a finite bound L' so that 

V qM-k n X kl ■ ■ ■ X k n < L> X kl--- X k n • 

k\---k n k\---k n 

The latter sum also contains a finite number of terms. 
Therefore a bound L" exists such that Y) x ki • ■ ■ Xk n < 

k\ • • • kn 

L" Y Xk- Consequently, there is a finite positive constant 
k 

L so that the mode-coupling polynomial is bounded by 
a linear one: 

Pq[Pe,x] < L £(x) , £(x)=5>, 

k 

(15a) 

for all q. Substituting this inequality into the 
TMCT map, Eq. (9b), one finds ]T)7^[P e ,x] < 

9 

M exp[— l/(L£(x))]. Elementary considerations show 
the following. There is a positive number £ max such 
that for all £ < £ max the inequality M exp[— 1/(L£)] < 
^exp(—1/L) is valid. This leads to 

^T g [P e ,x] < £(x) ex P (—l/i) , 0 < ^(x) < (max • 

(15b) 

The inequality (15b) is the basis for the proof of the 
theorem. The proof is done indirectly. Let us assume 
that there is a state P e with an associated fixed point 
f e = f(P e ) which tends continuously to zero for £ —> 0. 
In particular £ E = Y f q tends continuously to zero. Con- 

9 

sequently, one can find a state P e close to P c such that 
< £max- However, this result contradicts the inequality 
(15b) since exp(— 1/L) < 1 and J07^[P e ,f e ] = £ e , which 

9 

follows from the fixed-point equation ( 10 ). 
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III. THE FIRST SCALING LAW 

In this section the TMCT equations of motion shall be 
reformulated for a treatment of states P e with e tending 
to zero. This is done as originally performed for MCT 
(see Secs. 6.1.1 and 6.1.3 of Ref. [2]). The remaining 
discussions shall be restricted to the TMCT. 

In the limit of e —> 0 the states P e tend towards the 
glass-transition singularity P c . A first ansatz shall be 
made by treating (</> q (t) — f q ) as a small quantity of order 
yf\e\. The equations shall be reformulated up to terms 
of order |e|. Terms of higher order shall be indicated by 
the symbol O e , i.e. lim e _yo O e /e = 0. Let us write 

Mt) ~ fq = s q 9q(t) , (16a) 

s q = ln (V/ 9 C ) > (!6b) 

where the amplitude s q follows from Eq. (lib) for x = f c 
The ansatz implies 


dependence on control parameters enters via the smooth 
functions of e. 

D q [P E ] = s c q {P q [P £ ,f c ]-P q [P c ,f c ]} . (19) 

Further reformulations for K q (z) and E q (t) are needed 
below and they are presented in an appendix. 

Equation (A5) together with the formulas (A6) - (A8) 
for I q (t) provide the reformulated equation of motion for 
g q (t). How to solve such an inhomogeneous equation has 
been discussed in the last paragraph of Sec. IIC. The 
array X in Eq. (14a) can be identified with the array 
g (t) in Eq. (A5). According to the ansatz it is of leading 
first order in the small parameter It contributes to 

Eq. (14c) the term g e (t)a, where the factor X in Eq. (14c) 
is denoted by g e {t). According to Eq. (A6), the array I(t) 
consists of three terms of second order in y/fe[. Thus, the 
contribution X q in Eq. (14c) does not contribute to the 
leading-order results, to which the following discussions 
are restricted. Therefore Eq. (16a) takes the form 


3fci (*)•■■ 9k n it) = O e , n> 3 . (17a) 

For liquid states, the ansatz is non-trivial since restric¬ 
tions of g q (t) to order y'jij imply restrictions of the ad¬ 
missible time interval for small as well as for large t. A 
second ansatz concerns the regular term in Eq. (6). It 
will be assumed that it can be neglected relative to the 
mode-coupling kernel : 

[z + 'Wq\/[rhq{z)\ = O e . (17b) 

The two preceding formulas shall be justified a posteriori 
at the end of this section in connection with the definition 
of a scaling limit. 

The reformulations start by substitution of Eq. (16a) 
into Eqs. (4) and (5) in order to obtain the relaxation 
kernel as a polynomial in g (t) : 


s q m q (t) = 1 + w q lk9k{t)+ 
k 

+ E W S P gk(t)g P (t) + D q [P e ] + O e . (18a) 

kp 


The convention made in Eq. (16b) implies that the lead¬ 
ing order of the expansion of s q m q (t) is unity. The first 
order coefficients are related to the critical stability ma¬ 
trix: 

w$ = s c q dJ\/df k \P\ f c ]sl = [l/fq\A q k . (18b) 

The second order coefficients read: 

Kl = s c q d 2 Pq/df k df p [P c , i c ]s c k s c p /2 . (18c) 

The coefficients are to be evaluated at the critical point. 
Smooth variations with e can be absorbed in O e . The 


4>{t) = i c + hg e (t) , (20a) 

with the critical amplitudes introduced by 

h q = s q a q , q = 1, ■ • • , M . (20b) 

The preceding reasoning is correct if and only if the 
solubility condition Eq. (14b) is obeyed. The first contri¬ 
bution to ^P q a*I q (t) is a smooth function which vanishes 
for e = 0. It is called the separation parameter a(e): 

a(e) = Y, a * q fqDq[Pe] = Ce, C> 0. (21) 

q 

The second equality sign holds for generic paths, and if 
terms of order (y/RI) 4 are neglected. The second con¬ 
tribution yields A g £ (t) 2 with A = J2 q , k p a q^ q ,kp a k a p- 

Substitution of Eq. (A7) and eliminating W^ k by us¬ 
ing Eqs. (18b), (13a), and (13c) one gets the so-called 
exponent parameter 


A = E a *qfq W Sp a ^P + (V 2 ) E S VK ■ ( 22 ) 

q,kp q 

It differs from the exponent parameter within MCT, par¬ 
ticularly due to the presence of the second term. The 
third contribution is a product A 1 LT~ 1 [zg e (z) 2 ](t). Us¬ 
ing Eqs. (A8), (18b), (13a) and (13c) one gets A' = 

J2q a qfq^q,k WqJ a k a p = 1 - 

As a result, the TMCT equations of motion are re¬ 
duced to the scaling equation 


cr(e) + A g E (t) 2 - d t 



t')g e (t')dt' = 0 


(23) 
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This is the same non-linear integro-differential equation, 
which was studied in MCT. Only those implications shall 
be cited, which are needed for the discussion below. 

For cr = 0, Eq. (23) is solved by the critical correlator 
g c {t) = (f 0 /t) a . This holds provided the critical exponent 
a is related to the exponent parameter A according to: 

r(l - a) 2 /r(l - 2a) = A , 0 < a < 1/2 . (24) 

The ansatz (17b) has eliminated the effect of the tran¬ 
sient dynamics on the solutions. Therefore, the time scale 
to cannot be determined analytically. It has to be chosen 
such that the critical correlator matches the transient. 
Let 5 ± (t) denote the solution of Eq. (23) for a = 1 (glass) 
and for a = — 1 (liquid). The critical correlator shall be 
chosen as inital condition 

lim 5 ± (t)t a = 1 . (25a) 

t —>-0 

Then a solution of the scaling equation Eq. (23) is given 
by: 

g e {t) = c a g ± (t) , g e {z)=c a t a g ± (z) . (25b) 

Here rescaled times t and rescaled frequencies z are de¬ 
fined by 

t = t/t a , z = zt a . (25c) 

This holds provided the correlation scale c a and the time 
scale t a are chosen as 

Ca = \fW\, ta=to/\a\ s , S = l/2a . (25d) 

The preceding equations formulate the first scaling law. 
They describe the evolution of the dynamics in terms of 
control-parameter independent shape functions g±(i) or 
g ± (z). The sensitive changes, which occur if the states 
P e approach the glass-transition singularity are described 
by the variations of the two scales c CT and t a . 

It remains the task to justify the ansatz, Eq. (17b). In¬ 
deed, one gets [z + iv q \/rh q (z) = [z/t a +iv q ]/[rh q (z)t a \ = 
O ei because (1 /t a ) = {l/to)C\e\ • |cr(e)|’ 7 , with g = 
(<5 — 1) > 0. Behind the derivation of the first scaling 
law, there is the scaling limit for a > 0 : 

lim [(t> q {ito-)-f q ]/{h q c a ) = g ± {t ) , q = 1, ■ • • , M . 

e— »0,0<£ fixed 

(26) 

For e —> 0, there appears via the correlation scale c a 
the small parameter y/jej. This was anticipated by the 
ansatz, Eq. (16a). But, the time shifts towards infinity 
as tt a . The scaling law views the dynamics on scale t a . 

The shape function for the glass exhibits arrest: 
lim^^ g + {t) = l/\/l — A. The asymptotic law for the 
fixed point reads: 


fq ~ fq + hq , a —> 0 + . (27) 

The square-root variation is the signature of a fold bifur¬ 
cation. The divergency for A —> 1 signalizes the approach 
to a higher-order singularity. The 4> q (t)-versus-t curve 
describes the crossover from the critical decay to arrest. 

The shape function for the liquid and t —> oo exhibits 
power-law variation quantified by an exponent b: 

g~(t) ~ —Bt b , t —> oo . (28a) 

This is von Schwcidler’s law. The von Schweidler expo¬ 
nent b is related to the exponent parameter A : 

A = T(1 + &) 2 /r(l + 25) , 0 < 6 < 1 . (28b) 

The constant B cannot be calculated analytically. Sub¬ 
stitution of these results into Eq. (20a) casts von Schwei- 
cller’s law into the form: 

0 (f) = f c - h(t/o b 

(29a) 

t' a = i'o/M'* , t' 0 = to/S 1 / 6 ,7 = (1/2 a) + (1/2 b) . 

(29b) 

It describes power-law decay below f c . There appears a 
new time scale t' a with t' a lt a oc l/|cr| x / 2b 

IV. THE SECOND SCALING LAW 

In this section all considerations are restricted to liquid 
states close to the transition point: a < 0, P e —> P c . 
Validity of the second scaling law of MCT will be proven 
for TMCT. Again, only the crucial steps will be given. 
For more details see Secs. 6.2.1 and 6.2.2 of Ref. [2]. 

A. The equations of motion 

A second scaling limit shall be analyzed, again justified 
a posteriori by its success. It is based on the introduction 
of a cr-sensitive time scale r, which diverges more strongly 
than t a : lim e _>.o t a /r = 0. Rescaled times and rescaled 
frequencies are introduced: t = t/r, z = zt. The ansatz 
for the searched for array reads 0(f) = F(f, er) = 0(f) + 
cr tp(t) + O t . The scaling limit concerns the approach of 
(7 to zero for fixed f and z. The leading-order result, to 
which the discussions shall be restricted, reads 

0(f) = 0(f) , 0(z) = T0(z) . (30) 

The correlators are given by shape functions 0 which are 
independent of the control parameters. The change of 
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the dynamics due to changes of the state P £ is caused 
solely by changes of the scale r. Substitution of <p(t) 
from Eq. (30) into Eq. (5) for the mode-coupling kernels 
implies the scaling law: 

m (t) = m(i) = T\P C , 0(f)] . (31) 

Changes of the coefficients V^ i ... kn {P e ) in Eq. (4) due 
to changes of P £ would contribute only to the corrections 

The first step towards closed relations of the rescaled 
quantities is based on Eq. (6). It reads K q (z ) = 
— (1 /t)/{{z/t + iv q )/(Ttt q ) + m q (z)}. Here we used 
m (z) = rm(z), which follows from Eq. (31). For the 
scaling limit r —> oo, but t fixed, one gets 

r K q (z) = Kq(z) = - l/rh q (z ) . (32a) 

The second step deals with Eq.(7b). One gets from 
Eq. (7b) with Eq. (32a) 

E q (z) = T/[m q (z)z 2 ] = rEq(z) , (32b) 

which implies 

E q (t) = E q {t)=LT- 1 [l/[fh q {z)z 1 ]{t) . (32c) 

Hence, <j) q (t) = exp[— E q (t)] = exp[— E q (t)] = <f> q (i). This 
proves consistency of this procedure and it establishes 
the second scaling law. The scale r cannot be fixed after 
the parameters fl q and u q for the transient effects have 
been eliminated. 

The initial values of the shape functions, say, f q = 
4 > q (t —> 0), are implications of the preceding equations, 
which do not depend on the scale r. One gets using 
Eqs. (7b), (32b) and (31): 

fq z) = exp [-Eq(t -A 0)] = exp[lim?-^ [zE q (5)]] = 
exp[lim 5 _ Voo [l/(m g (5)i)]] = exp[-l/m g (f -A 0)] = 
exp[— l/P q [P c , fM]]. Equations (9b) and (10) imply that 
fW is a fixed point for P = P c , Generically, one gets 

4>(t -A 0) = f c . (33) 

The second scaling law does not deal with the decay of 
correlations below unity, but it deals with the decay be¬ 
low r q . 

B. von Schweidler’s law 

Equation (33) suggests to solve the equations of motion 
for the limit t tending to zero. Writing 

(34) 


the solution can be obtained using g q (t) as small quanti¬ 
ties. Substituting these expressions into Eq. (31) yields 
the expansion up to terms in second order: 

s c q rh q (t) = 1 + ^2 W q l l9k (t) + J2 W q 2 k P 9k C t)g P ( t) . 

k kp 

( 35 ) 

The coefficients W q ^ an d W q 2 ^ p are given in Eqs. (18b) 
and (18c). This expansion in terms of g q (t) is identical to 
that given in Eq. (18a). Here, the procedure is simplified 
since the term D q [P e ] does not occur. The further re¬ 
formulations of the equations of motion can be obtained 
following the ones explained in the appendix. One ends 
up with Eq. (A5) with the inhomogeneity from Eq. (A6), 
but with D q [P e ) = 0. 

The further discussion follows that presented in 
Eqs. (20a), (20b), (21) and (22) for the derivation 

of the scaling equation (23). One finds in leading order 
g(f) = a g(t), where a is the right eigenvector of the crit¬ 
ical stability matrix specified in Eqs. (13a)- (13c). Func¬ 
tion g(t) is a solution of Eq. (23) specialized for er(e) = 0. 
A solution consistent with the requested short time be¬ 
havior reads g(t) = —t b . The positive exponent b is de¬ 
termined by Eq. (28b), where the exponent parameter A 
is noted in Eq. (22). For small rescaled times t = (t/r), 
there holds von Schweidler’s law: = f c — ht 6 . Up 

to here the scale r is not fixed. It is determined by the 
matching request. If, and only if, the long-time asymp¬ 
tote on scale t„ is identical with the small-time asymptote 
on scale r the second scaling law is a continuation of the 
first one. Equation (29a) requires that r = t' a . 


V. A SCHEMATIC MODEL 

Some features of the results shall be demonstrated by 
considering a model, which is simplified with the mere 
intention to reduce the numerical effort. First, the array 
of M components are chosen for M = 1, i.e. , only a 
single function <j>(t) is used. Second, the mode-coupling 
polynomial is chosen as 

m[P, x\ = v\x + V3X 3 . (36) 

The two coupling coefficients v\ and V 3 are used as con¬ 
trol parameters. A state is presented as a point of the 
first quadrant of the v± — v$ plane: P = (iq, V 3 ), v± > 0, 
V 3 > 0. This mode-coupling polynomial is a schematic 
version of the corresponding polynomial for a spin glass 
transition [16]. 

Within TMCT the fixed-point equation (10) reads v\ + 
v 3 / 2 = l/[/hi(l//)]. The transition singularities are 
given by the request that the critical stability “matrix” 
from Eq. (12b) is unity: v^ + 3v%f c2 = l/[/ c (ln(l// c )) 2 ]. 
These two equations imply formulas for the critical point 
P c = (t^Ug) as function of the arrested parts f c : 


= fq+S q g q {t) 



Vi 

FIG. 1. (Colour online) Glass-transition diagram for a 
schematic model specified by a mode-coupling function 
J-[P, x ] = V1X+V3X 3 within MCT [19] and TMCT. The dashed 
line marks the positions of continuous and the solid lines rep¬ 
resent positions of discontinuous transitions. Both dots mark 
the state point (vi,V 3 ) with maximum exponent parameter 
A ma x = 2/3 for TMCT (full dot) and A ma x = 1 for MCT 
(open dot). The solid line for MCT terminates at the open 
dot. 


v c 3 = [3Gi(/ C ) - G 2 (/ c )]/2, u 3 c = [G 2 (/ c ) - G 1 (/ c )]/[2(r) 

(37a) 

Gi(/ C ) = l/[m G 2 (/ c ) = l/[/ c (s c ) 2 ], s c = ln(l// c ) 

(37b) 

The exponent parameter from Eq. (22) becomes 
A = s c [2 - (3/2)s c ] . 

(37c) 

The glass transition diagram calculated within TMCT 
is presented in Figure 1. The corresponding diagram cal¬ 
culated within MCT [19] is added for reasons of com¬ 
parison. The figure demonstrates that the trend to glass 
formation is stronger in MCT than in TMCT, as already 
noted in Ref. [15]. For = 0, the TMCT value for w 3 ex¬ 
ceeds that for the MCT by about 20%. For wj = 0, both 
theories yield A = 1/2, i.e. the von Schweidler exponent 
b equals unity. For this case, there is no stretching for 
the decay which is described by the second-scaling-law 
process. 

Increasing vf, the value u 3 decreases. Within MCT, 
this decrease is connected with an increase of A. The ex¬ 
ponent parameter reaches unity at a higher-order glass- 
transition point(open dot in Figure 1). This increase is 
connected with a decrease of the exponent a and b to¬ 
wards zero. The stretching increases up to the point, 
where the range of validity of both scaling laws shrinks 
to zero. Within MCT, also a transition line(dashed line 


in Figure 1) occurs, which is connected with a continuous 
variation of the arrested parts. This line terminates at a 
crossing point with the above mentioned transition line of 
discontinuous changes of the arrested part. This crossing 
implies the existence of glass-to-glass transitions. 

The TMCT does not exhibit any of the subtle phenom¬ 
ena described in the preceding paragraph. The w 3 -versus- 
v\ curve decreases monotonously through the quadrant 
if f c decreases from e -1 / 3 (at = 0) to e -1 (at vf = e). 
With increasing the exponent parameter increases up 
to a maximum value A max = 2/3 (at v\ = 9e 2//3 /8) . 
This maximum value is connected with f c = e~ 2 / 3 . For 
increasing u 3 further, A decreases. Consequently, there 
exists a minimal value of the von Schweidler exponent 
bmin ~ 0.70 corresponding to A max = 2/3. 


VI. SUMMARY AND CONCLUSIONS 

We have elaborated on similarities and differences be¬ 
tween the novel mode-coupling theory (TMCT), recently 
been proposed by Tokuyama [13], and the conventional 
mode-coupling theory (MCT). Our focus has been on 
the glass-transition singularities and the validity within 
TMCT of both scaling laws found within MCT. 

A very important difference between both theories 
is the absence of a continuous glass transition within 
TMCT with mode-coupling polynomial (4). This is re- 
21 lated to the fact that all the derivatives dT q [P, x]/dxk 
vanish at x = 0 within TMCT, in contrast to MCT. 

As an illustration the phase diagram of a schematic 
model has been determined analytically. Figure 1 shows 
that there is a transition line within MCT at which con¬ 
tinuous glass transitions occur, in contrast to TMCT. 
This transition line terminates at the MCT-transition 
line of discontinuous transitions. Between the endpoint 
(open dot in Figure 1) of the latter transition line and 
the crossing point there occur glass-to-glass transitions. 
These phenomena do not exist for TMCT. Another im¬ 
portant difference is the absence of higher order glass- 
transition singularities for the schematic model within 
TMCT, because the exponent parameter is always be¬ 
low unity. These results challenge within TMCT in¬ 
vestigations of microscopic models like short-range at¬ 
tractive colloids. For such liquids MCT predicts higher 
order glass-transition singularities, glass-to-glass transi¬ 
tions and re-entrant behavior [20-22]. 

The first scaling law found for MCT implies for the 
relaxation of the correlators two time power laws with 
diverging time scales. These involve exponents which 
are determined by the exponent parameter A. We have 
proven that TMCT leads to the same scaling equation 
but with a different expression for A. Consequently, the 
long-time relaxation obtained within TMCT also exhibits 
two power laws with diverging time scales, however, with 
different numerical values for the corresponding expo¬ 
nents. Explicit formulas for the evaluation of the ex¬ 
ponent parameter and for the amplitudes of the singular 
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parts of the correlators have been obtained. The second 
scaling law within MCT proves the existence of the su¬ 
perposition principle, e.g. for the time and temperature 
dependence of a correlator. This scaling law also holds 
within TMCT. 

The kind of scaling laws originally found within MCT 
were not known before in the field of ’’Nonlinear dynam¬ 
ics” , where bifurcations of dynamical states are inves¬ 
tigated. The validity of both scaling laws within MCT 
and TMCT proves some robustness of the dynamical fea¬ 
tures of the bifurcation scenario which describes the ideal 
liquid-to-glass transition. 
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Appendix A: Reformulation of the TMCT equations 
of motion 

Laplace transformation of Eq. (18a) leads to 

~s c q zrh q (z) = 1 ~zJ2 W^g k (z) + D q [P e ] 
k 

~z E W^ p LT[g k (t)g p (t)\(z) + O e . (Al) 

kp 

The ansatz in Eq. (17b) means that those terms [z+iv q \ 
in Eq. (6) which specify the transient dynamics can be 
dropped. Substituting Eq. (Al) into Eq. (6) yields 

k q {z)/{s c q z) = 1 + zJ2w^g k (z) - D q [P e ] 
k 

+z E w q 2 kp LT lffk(t)g P (t)](z) 

kp 

+ {zJ2 W S9k(z)} 2 + O e . (A2) 

k 

Using Eq. (A2) one obtains a corresponding expansion 
of E q (z) = — K q {z)/z 2 (cf. Eq. (7b)). Then transforma¬ 
tion into the time domain leads to 


Substituting this expression into the formula (f> q (t) = 
exp {—E q {t)) (cf. Eq. (7b)) and expanding the exponen¬ 
tial function, one arrives at 

<t>S)/f q = 1 + s c q E W$g k (t) + a° q D q [P e ] 

k 

+S q E W q%9k{t)g p {t) 

kp 

S q E Klw$LT- l [zg k {z)g p {zm 

kp 

+ \(K) 2 E W S W ^ 9 k(t)g p (t) + O e . (A4) 

kp 

Equation (16a) reads <t> q {t)/f q - 1 = [s c q J f q \g P {t). As a 
result, one obtains the desired reformulation of the equa¬ 
tions of motion 

E^ “ A q,k\9k(t) = I q {t) , (A5) 

k 

with (Ay k ) the stability matrix (Eq. (12b)) at the critical 
point. 

The M components of the array 1(f) are 

IS) = f q D q [Pe\ + E ^ q ,k P 9k(t)g P (t) 

kp 

+ E N^ kp LT- l [zg k (z)g p (z)} (f) + (A6) 

kp 

with 

= ’ ( A7 ) 

N c q , kp = f c q W$W$ , (A8) 

where and W^ 2 kp are given by Eqs. (18b) and (18c), 
respectively. 


E(t)/S c q = 1 - E <iW) - D q [P e ] 

k 

- E W q 2) kp9kit)g p {t) 

kp 

- E LT ~ X i-9k{z)g P {z)] (t) + O s . (A3) 

kp 
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